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This is a write-up of a poster presented in the international conference Quark Matter 2004 at 
Oakland, CA. We will report centrality dependence of chemical freeze-out temperature (T c h), light 
quark chemical potential (fi q ), strange quark chemical potential (fi s ), and strangeness saturation 
factor (7 S ) in Au+Au collisions at -^/sjvjv = 130 and 200 GeV. A systematic study for combination 
of ratios for chemical freeze-out fit is studied and we have found small combination dependences. 
The results show j a increasing with centrality but the other parameters have less sensitivity to the 
centralities. 
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I. INTRODUCTION 

The statistical model approach is established by analy- 
sis of particle ratios of the high energy heavy ion collisions 
in GSI-SIS to CERN-SPS energy [iBllREllll and el- 
ementary collisions (e + e~, pp, and pp) [alii- The model 
describes the particle ratios by the chemical freeze-out 
temperature (T c h), the chemical potential (/x), and the 
strangeness saturation factor (7^ ) . Many feature of the 
data imply that a large degree of the chemical equili- 
bration may be reached both at AGS and SPS energies 
excepting strangeness hadrons. There are the four most 
important results. 

1. At high energy collisions the chemical freeze-out 
(inelastic collisions cease) occurs at about 150-170 
MeV and it is 'universal' to both elementary and 
the heavy ion collisions; 

2. The strangeness is not fully equilibrated because 7^ 
is 0.5-0.8 [a, EEl (if strangeness is in equilibration, 

7.=i 0); 

3. The kinetic freeze-out (elastic scatterings cease) oc- 
curs at a lower temperature ~ 100-120 MeV; 

4. The compilation of the freeze-out parameters [Tol ] in 
the heavy ion collisions in the energy range from 1 
- 200 A- GeV shows that a constant energy per par- 
ticle (E) I (N) ~ 1 GeV can reproduce the behavior 
in the temperature-potential (T h — fi B ) plane |10| . 

We have many hadron yields and ratios including multi- 
strange hadrons as a function of centrality in Au+Au 
collisions at ^/snn = 130 and 200 GeV. They allow us 
to study centrality dependence of chemical freeze-out at 
RHIC energy. 

II. MODEL 

We adopt the grand canonical ensemble for a chemical 
freeze-out model based on Ref. [Ill and the model used in 



Ref. Gnami. We assume no difference between u and 

d quark chemical potentials. Therefore, the hadron gas 
is described by a chemical freeze-out temperature (T c h), 
light (u and d) quark potential (n q ), strange quark po- 
tential (/is), and strangeness saturation factor (7s )• The 
density of a particle i in the hadron gas is given by: 

Pi =7 S <5+S_> « J^ ch 3 (jp-Y K 2 ( mi /T ch ) \ q Ql V 

(1) 

where m$ is the mass of the hadron i, gi is the number 
of spin-isospin degree-of-freedom, Ki is the second-order 
modified Bessel function and, 

X q = exp(fi q /T ch ), X s = exp((j, s /T ch ). 

The potential fi q is for u/d/u/d quarks, and fj, s is for 
s/s quarks. The /i q is a third of baryon chemical po- 
tential \lb. Qi ans correspond to the net number of 
valence u/d quarks (Qi — (u— u + d — <i) i ), and s quark 
(sj = (s — s).) of particle species i, respectively. The 
RHIC data we used are from mid-rapidity region, there- 
fore we does not apply any conservation low and treat 
all of chemical freeze-out parameters as free parameters 
for the fit. The factor 7 S (0<7 S ) is introduced to take 
account of possible incomplete chemical equilibration for 
strange particles 1]. 7 S =1 signifies full strangeness chem- 
ical equilibrium and 7 S >1 shows over strangeness sat- 
uration. The power factor of 7 S is total number of s 
and s quarks in the particle. It is not clear whether 7^ 
should be unity or not. On the other hand, calculations 
in Refs. |a,[l3>[il,ll9 support 7 S <1. Therefore, we adopt 
7 S as a free parameter. Note that we use the Boltzmann 
approximation for all hadrons except pions, where the 
Bose distribution is applied. The particle densities are 
computed by Eq. for the hadron gas including known 
resonances up to mass of 1.7 GeV/c 2 , and the effect of 
broad resonance mass widths is taken into account |llj| . 
The resonance decay to lower mass hadrons after chemi- 
cal freeze-out is also taken into account. 

There are technical caveats in above approaches. The 
application of some statistical model requires that the 
measurement is done over the whole phase-space, i.e., 
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with 4-7T particle yields. This is because conservation laws 
that apply to the collisions are only valid for global mea- 
surements, not locally. In all experiments the coverage 
of the phase-space is limited. Therefore, extracting the 
yields measured within limited phase-space to Anx yields is 
almost impossible, especially for the collider experiments. 
In addition, we should point out that the non-local effect 
is common to all global conservation laws (baryon num- 
ber, charge, and so on), not only to the strangeness. 

We adopt the static model that assumes full phase 
space to limited phase space, that is, limited rapidity 
region. Ref. [l7j has studied effects of excluded vol- 
ume correction, limited rapidity range, transverse flow, 
and longitudinal expansion to the ratios. Various effects 
which can severely distort the momentum spectra of the 
particles produced cancel out in such ratios ^?|]- This 
addresses that the static model can be applied to the 
particle ratios in the limited rapidity range. 



III. ANALYSIS 



Tabled lists the centrality definition and references for 
the value of dN/dy and the ratios at mid-rapidity for 
each experiment. The centralities are defined by the in- 
dependent detector and different acceptances among the 
experiments. The common parameter of the collision cen- 
trality is an average of number of participant ((N part }). 
Therefore, we use (N part ) to discuss the centrality depen- 
dence of the chemical freeze-out parameters. 

We select two kind of centrality binning from STAR 
data. The first one is used for data set from charged pion, 
kaon and proton (centrality selection used in Ref. [2^.l29| 
for 130 GeV and Ref. :38] for 200 GeV). The second 
one has smaller number of bins but more particle ra- 
tios (centrality selection used in Ref. [2^| for 130 GeV 
and Ref. ^(j f° r 200 GeV). The other experiments have 
different centrality bins, therefore we need to estimate 
particle ratio for the centrality bins selected. At first, we 
interpolate the dN/dy of each particle assuming linear 
scaling of (N part ) between two data points. After we ob- 
tained dN/dy values for the (N part }, the particle ratios 
is computed from the dN/dy values. 

The chemical freeze-out model fit is applied to the par- 
ticle ratios. A program package MINUIT |43j is used for 
the fit to find \ 2 minimum points and the error of fit pa- 
rameters is estimated as y 2 minimum + 1. The chemical 
freeze-out parameters seem to be sensitive to combina- 
tion of particle ratios as discussed in Ref. Q . Hence we 
check the following six combinations of particle ratios for 
the fit: (1) tt, K, and p; (2) tt, K, p and A; (3) tt, K, p, 
A, 0, and H; (4) tt, K, p, A, K* , (j>, and H; (5) tt, K, p, 
A, <fi, S, and J7; (6) tt, K, p, A, K* , <fi, S, and fi. 



TABLE I: List of possible dN/dy and ratios from the experi- 
mental data in the Au+Au collisions at y / sjvjv = 130 GeV and 
200 GeV. The parenthetic particle name ((A* ) and (K^)) 
means average of the particle and the anti-particle. The se- 
lection of event fraction is defined by (I) the multiplicity of 
charged particles from the Multiplicity Array [l8|. (II) the 
energy deposit on Zero Degree Calorimeter (ZDC) and the 
number of charged particle measured by Beam-Beam Counter 
(BBC) [ Tg || . (Ill) the signal of the scintillator paddle coun- 
ters |20| . (IV) the energy deposit on ZDC and the signal of 
Central Trigger Barrel (CTB) (V) the number of recon- 
structed charged tracks (77 < |0.75|) |22| . (VI) the number of 
reconstructed charged tracks (77 < |0.5j). Additionally, some 
references show the number of participant and the dN/r/ of 
negative charged hadrons for the event fraction. We marked 
them in the list by {N par t) and dN h - /dr). 
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IV. RESULTS AND DISCUSSION 

Figure 2] shows centrality dependence of chemical 
freeze-out temperature T c h, light-quark potential fi q (= 
/xs/3 where /xs is Baryonic chemical potential), strange 
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FIG. 1: Chemical freeze-out temperature T c h, light-quark 
potential /j, q , strange quark potential fi a , and strangeness sat- 
uration factor 7 S as a function of average number of partic- 
ipant. The different markers show different combination of 
particle ratio used for the fit. 
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FIG. 2: Comparison of fit results to the data in 130 GeV 
Au+Au central collisions ((iVp ar t)=317). All of data from 
dataset (6) shown in the plot are used for the fit. The point 
with marker is experimental data. In the top plot, the model 
calculations of ratio from the chemical freeze-out parameters 
obtained are shown horizontal line. The variation in the error 
of the parameters is shown as shaded box. The bottom plot 
shows a deviation of data to the model, (r exp — r mo dei)/ Ar exp , 
where r exp is ratio from data, r m odei is ratio by model calcu- 



quark potential fj, s , and strangeness saturation factor 7 S 
as a function of (N part ) in 130 and 200 GeV Au+Au col- 
lisions. The markers correspond to combination of ratios 
used for fit as shown in the figure. Table [H] lists all of 
values shown in Fig. ^ 

The chemical freeze-out parameters vary with combi- 
nation of particle ratios and shows larger variation in 130 
GeV Au+Au than 200 GeV. Taking care the variation, 
T c h, (J-q, and /i s seems to be constant with centrality in 
both 130 GeV and 200 GeV Au+Au collisions. The 7 S 
is slightly increasing with centrality in 130 GeV when 
we focus the results by the fit to ratio from it, K, p, A, 
K* , </>, 5 (,and fi). On the other hand, we can see the 
rise of 7 S as a function of (N part ) in any combination 
of ratios in 200 GeV Au+Au. The j s is about 0.5-0.7 
in the heavy ion collisions at SPS energies ^(|, e + e~ at 
LEP, and pp at SppS 0, 13 • Only Au+Au central colli- 
sions at RHIC energy (^/snn = 130 and 200 GeV) reach 
fully strangeness equilibration. The rise of 7 S seems to 
be rapid until around (N part ) = 100-150 and then slowly 
increasing. 

The T c h, fi>q, fi>s seem to be independent as a function 
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FIG. 3: Comparison of fit results to the data in 200 GeV 
Au+Au central collisions ({N par t) =322). All of data from 
dataset (6) shown in the plot are used for the fit. The legend 
is the same with Fig. [5] 



from combination dependence, T c h is 150—170 MeV, fx q 
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is 8—11 MeV, and (x a is -3—3 MeV. Small but finite \x q 
reflects net baryon density is small but not zero around 
mid-rapidity. The /i s is consistent with zero. It suggests 
that the chemical freeze-out occurs at phase boundary of 
a hadron gas. 0, US • 

The chemical freeze-out parameter varies by data com- 
bination. In the central 130 GeV Au+Au data, T c h, 
fx q , fi s , and 7 S change about 10 MeV, 4MeV, lMeV and 
0.1, respectively. The deviation is small and consistent 
within error for all of parameters. The variation becomes 
smaller in the central 200 GeV Au+Au data. It is about 
5 MeV, 2MeV, 2MeV and 0.06 for T ch , n q , fi s , and j s , re- 
spectively. The results from both energy shows that the 
7 S increases when we adopt multi-strangeness in the fit at 
a centrality bin. The "f s is sensitive to multi-strangeness, 
therefore the constrain changes stronger with increasing 
number of strangeness data point in the fit. However, it 
is not clear how the number of strangeness play the role 
in the constrain. We need to check this point with more 
data that will come soon from RIIIC run4. 

Figure[3and|nishow a comparison of data to the model 
calculation. The data is included multi-strangeness and 
it is all of possible data sample at this time (dataset (6)). 
The model calculation shows good agreement in both en- 
ergy, excepting K*. In the 200GeV, the (A + A)/7r~ value 
from model also differ to the data. Note that the data 
A + A is obtained from a slide of the conference presen- 



tation from Ref. |4lJ and the value seems to be changed 
in Quark Matter 2004 conference, therefore the results 
might be changed and we need to check in future. 



V. SUMMARY 

We study centrality dependence of chemical freeze-out 
parameters in ^snn = 130 and 200 GeV Au+Au colli- 
sions. The T c h, Mgi Ms show independence as a function 
of (Np ar t). The most remarkable point is about central- 
ity dependence of j s . The result shows 7 S increasing 
with centrality and reach fully strangeness equilibration 
in only Au+Au central collisions at RHIC energy. The 
rise of j s seems to be rapid until around (N part } = 100- 
150 and then slowly increasing. The ratio combination 
dependence is also studied in both energy and we have 
found small combination dependences. The T c h, /J. q , /i s 
seem to be independent as a function of (N part ). On 
the other hand, the 7 S increases when we adopt multi- 
strangeness in the fit at a centrality bin. 
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TABLE II: Chemical freeze-out parameters for each average number of participant ((N P art))- The data set is related to 
combination of ratios used for the fit: (1) n, K, and p; (2) tt, K, p and A; (3) ir, K, p, A, <j>, and E; (4) tt, K, p, A, K* , (j>, and 
3; (5) n, K, p, A, 4>, 3, and fi; (6) n, K, p, A, K* , 4>, 3, and fi. The parameters are obtained by fit of chemical freeze-out 
model to the data. The error of parameters are defined as x 2 minimum + 1. 
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(31 


317 


168 fi+fi 1 

J-Uu. vj i vj . ± 


13.2+1.4 


2.0+1.5 


980+0 059 

VJ • i J O V J 1 VJ . VJ • J o 


11.5/11 




211 


1 65 fi+5 9 


11.9+1.8 


1.7+2.5 


930+0 060 

VJ • i J » J V J 1 V J .VJVJVJ 


8.4/11 




64 


1 65 5+5 6 

lUtJ. 'J 1 'J • VJ 


11.5+1.8 


0.8+2.5 


850+0 055 

VJ • LJVJVJ 1 VJ . VJ 'J 'J 


17.2/11 


(51 


317 


169 2+fi 

1U J. ZJ_I_U. VJ 


13.2+1.4 


1.9+1.5 


984+0 057 

VJ • Jul 1 VJ . VJ < J ( 


11.7/13 


(61 


317 


169 1+6 1 

_1_ VJ tj • J 1 — Ui 1 


13.2+1.4 


2.0+1.5 


972+0 057 

Uti/ 1 ^ — 1 — *J • w 'J I 


13.9/14 


(11 


351 


152 4+3 3 


7.9+2.1 


1.5+4.6 


963+0 063 

VJ • i J VJ tj 1 VJ . VJ VJ (J 


1.2/ 3 




298 


152 1+3 3 

_L tJZj m A 1 \J m *_J 


7.7+2.0 


2.1+4.6 


942+0 062 

VJ • * / 1 1 V / . \ J V ; Z- 


1.1/ 3 




234 


151 fi+3 2 

J_tJ _L . V J 1 (J . Zj 


7.4+2.0 


0.3+4.6 


927+0 060 

VJ . iJZj 1 1 VJ.VJVJD 


1.4/ 3 




166 


151 2+3 2 


7.6+2.0 


0.6+4.5 


902+0 059 

VJ . iJVJZj_I_VJ . VJ < J ! J 


1.9 / 3 




115 


151 2+3 2 

_L xj _L ■ ZJ 1 ■ ) - —t 


7.3+2.0 


0.6+4.5 


884+0 058 

VJ . VjLJti—I—VJ . VJ < J O 


2.3/ 3 




76 


150 3+3 1 

_l_ xj vj ■ ' / — l — *J • 1 


7.2+2.0 


-0.5+4.5 


853+0 056 

U • "J — 1 — VJ • \ i ' 1 \ t 


3.1/ 3 




47 


149 5+3 1 

1t:J. 'J 1 (J . _L 


6.8+2.0 


-0.6+4.5 


822+0 054 

VJ . LJZjZj_I_VJ . VJ ' J ^ 


4.9 / 3 




27 


147.4+2.9 


6.4+2.0 


-1.3+4.5 


776+0 051 

VJ . 1 1 VJ 1 VJ . VJ'J X 


7.2/ 3 




14 


143.1+2.7 


6.5+2.0 


3.1+4.5 


702+0 047 

\J . 1 VJZj_I_ VJ . V7t: 1 


11.4/ 3 


(21 


322 


156 5+2 9 

1 T J \ J • ' ' — 1 — Z. ■ . ■ } 


9.3+1.2 


3.0+2.6 


999+0 056 

U ■ ij ij t/ — I — vj ■ u rj vj 


10.3/ 4 




234 


156 5+3 5 

J-VJVJ . <J 1 'J • Xj 


9.6+2.4 


0.0+4.7 


987+0 059 

VJ . JU 1 1 VJ . 1 J < J i J 


13.1/ 4 




140 


156 3+3 6 

J-VJVJ . tj 1 tj . VJ 


9.7+2.4 


0.3+4.7 


952+0 057 

KJ • ij <J _l_ VJ . VJ tj 1 


15.1/ 4 




76 


155 3+3 5 


9.5+2.4 


-0.9+4.7 


910+0 055 

VJ. ij J. VJ 1 VJ . V J < J < J 


16.5/ 4 




29 


152 4+3 3 

J_tJZj .1 1 'J . ' 7 


8.7+2.3 


-1.8+4.7 


833+0 050 

KJ • VjfJfJ_l_VJ . V J < J VJ 


20.6/ 4 


(31 


322 


1 57 8+2 7 

J-U 1 . U 1 ZJ . 1 


9.7+1.2 


3.5+2.4 


1 01 1+0 044 

_L ■ VJ J_ J 1 VJ ■ U ± ± 


12.1/ 7 




234 


158.5+3.2 


10.9+2.2 


2.3+3.8 


0.990+0.046 


14.4/ 7 




140 


158.6+3.2 


11.4+2.2 


2.6+3.6 


0.968+0.045 


16.9/ 7 




76 


157.2+3.0 


10.8+2.1 


0.9+3.6 


0.920+0.042 


17.4/ 7 




29 


153.4+2.9 


9.3+2.1 


-1.5+3.8 


0.850+0.043 


21.1/ 7 


(4) 


322 


156.3+2.7 


9.4+1.2 


3.4+2.4 


1.020+0.045 


19.5/ 8 




234 


157.2+3.1 


10.3+2.1 


2.0+3.7 


0.999+0.046 


19.3/ 8 




140 


157.8+3.0 


10.9+2.1 


2.5+3.6 


0.965+0.044 


18.6/ 8 




76 


156.8+3.0 


10.6+2.1 


0.9+3.6 


0.922+0.042 


18.3/ 8 




29 


153.2+2.8 


9.2+2.1 


-1.6+3.8 


0.852+0.043 


21.4/ 8 


(5) 


322 


157.9+2.7 


9.7+1.2 


3.2+2.3 


1.016+0.042 


12.5/ 9 


(6) 


322 


156.5+2.6 


9.4+1.2 


3.1+2.3 


1.027+0.043 


19.9/10 



130 GeV Au+Au 



200 GeV Au+Au 
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